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Abstract — In this paper, the most specific problems related to 
solar radiation and its conversion are treated. A numerical 
study of the dynamic operation of a solar water heater is 
performed. We focused on the choice of the weather data in a 
given location with taking an interest in the theory. A MATLAB 
program was designed to model the system’s operation. 

For this purpose, our approach consists in establishing the 
energy balance on each element of a solar water heater in 
dynamic operation. Then, we identify the weak points of the 
system and then looking for the parameters involved for a better 
design of the system adopted. Based on the different graphs, a 
good design of a solar water heater in a given location is made 
and this is done by selecting all critical parameters resulting in 
the highest efficiency or even when its variation is negligible so 
that the use of this parameter can save for the overall system. 

We gave a theoretical analysis that can lead to tangible results 
by following a computer work and essentially by developing a 
numerical code able to calculate the solar radiation in a given 
location based on the values found for calculating efficiency, the 
storage power and other adopted parameters. 

To improve the dynamic operation of the system as well as the 
economic aspect, the influence of different parameters on the 
efficiency of the overall system such as: the critical surface, the 
ratio of the mass of water to the surface of panel, the choice of 
flow of coolant, the coil length in the heat storage tank and the 
choice of coolant is investigated. 

Index Terms — solar thermal panel, solar water heater, solar 
energy 
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Ta 

Ambient 

temperature 

U 

Conductance of sensor 
[W/m 2 K] 

Tint 

Temperature around 
the balloon 

M 

Water mass into the 
balloon [kg] 

T u 

Usable temperature 

Gn 

Global lighting [W] 

To 

Networks 

temperature 

L 

Usage rate [kg/s] 

T f 

Temperature at the 
input of the panel 

A 

Area of panel 

Te 

Temperature at the 
output of the 

balloon 

Rt 

Thermal resistance of 
the tube [k/W] 

Ti 

Temperature at the 
input of panel 

Rb 

Thermal resistance of 
the coil [k/W] 

T s 

Temperature at the 
output of the 

balloon 

Qb 

Power stored in the 
balloon [W] 

T b 

Temperature of the 
storage tank 

Qu 

Useful power 

recovered by the 
balloon [W] 

F 

Correction 

coefficient 

Qp 

Total lost power [W] 

Fr 

absorber 

Qc 

Consumed power [W] 



Conductance factor 



F 

Efficiency of fin 

a 0 

Coef. absorption of 
the absorber 

(UA) b 

Conductance 
coefficient of the 
balloon [W/K] 

U 

coefficient of 

transmission of the 
glazing 


i. Introduction 

Solar energy has been used for a long time to produce 
heat. Although the expertise in this sector is considerable, 
nowadays an important renewal on the theoretical and 
experimental levels is known with the current 
technological development. Also, many studies have been 
conducted on the solar sensors by Sefeir et al [1] and 
Bernard [2]. 

For a long time, the specialists of solar thermal energy are 
interested by the modeling because it allows the performance 
of the facilities prediction. The numerical studies are fast, 
relatively cheap and can, in some circumstances, replace the 
experimental studies which are long and expensive. The 
majority of simulation softwares were developed to achieve 
the productivity of a solar water heater. 

The main objective of this work is to develop a numerical 
code able to calculate the solar radiation at a given location. 
By using the values found, we can also calculate the 
efficiency, the storage power and other settings adopted for 
the best choice of elements of the system. 

II. SIMULATION OF A SOLAR HOT WATER 

The modeled system consists of a glazed flat-plate panel 
supplying heat to a storage tank (Figure 1). 



Figure 1 : Scheme of a system for producing hot water 


Model adopted for the solar collector 

In this study, we used a glazed flat plate sensor (Figure 2) and 
an absorber finned tubes (Figure 3). The absorber has a length 
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L a , a width l a , a thickness e, the outer diameter of the tube is D e 
and “w” is the distance separating two successive tubes. 



Figure 2 : Glazed flat panel 



o,- l-- 

Figure 3 : Absorber finned tubes 


At the level of a portion of length L 2 between the 
output of the balloon and the input of panel: 

r \ 

(Ti — T int) = (Ts — T int)exp 


L2 


Rt.rn* .Cp y 
At the level of length coil L b in the balloon: 

Lb 


(:Ts - Tb ) = (re - 77?) exp 


Rb.rn .Cp 


- At the level of the storage balloon: 

In such a system, the water in the storage tank will be 
generally laminated and the hot water is located in the upper 
part. 

Therefore, the heat balance becomes: Qu = Qb + Qp + Qc 

M.C dT b r / \ -i 

ft = j - =A.Fr.R [ra 0 .Gn-U{T b -T a )]-Q p -Q e [2] 
at 


The amount re p resents power transferred to 

dt 

the continuous fluid into the balloon. R with the control 
function such as R = 1 when 

Tf > T t ; R = 0 when 7} < T r 

The lost power is expressed by: 

Q p =(U.A) b (T b -Tj 


III. HEAT BALANCE SYSTEM 

A. At the sensor level 


The consumed power is expressed by: 

Qc = L.C r (T„-T 0 ) 


The Useful power Qu collected by the panels is given by the 
equations of HOTEL-WHILLER and BLISS [3]. 

Qu = A.Fr.[ra 0 .G n -U(Ti-Ta)] 1 

Qu = m'.Cp. (Tf — Ti ) j 


Avec: Fr = 


F' = 


GmCp 


U 


1-exp 

1 

U 


U.F' \ 
Gm.Cp J 


W 


1 


F = 


th 


U[(w-De).F + De\ 
(W-De 


+ l.Raf 


my- 


w — De 
2 



B. At the pipe level 


- At the level of a portion of length Li between the 
output of the sensor and the input of the balloon: 

We isolate a portion of tube length dx and we establish the 
balance of the fluid. Between the times t and t-vdt , after 
integration, we obtain: 


(Te — T int) — (Tf — T int )exp 


LI 


Rt.rn .Cp 


As we can see, there are some losses corresponding to the 
fluid flow and its contact with the walls of the pipe. 


IV. SOLVING EQUATIONS 

In fact, the parameters Ti, Ti, Te, Tb, Ts, Fr, a, Gn, Ta et Tint 
depend on the time. In our study, F r and a are assumed 
constant, although these two parameters depend on several 
parameters including: the nature of the glazing, the nature of 
the absorbing plate, the speed of wind and the external 
temperature. The value of Fr was taken equal to 0.8 and that 
of a= 0.46. 

To solve equations [2] which governing the dynamic 
operation of our system, the method of EULER is used. This 
method is characterized by the use of finite differences to give 
us sufficient accuracy. The integration of these equations is 
usually done with the time step on the order of the hour. 
Therefore: 

U + i)=U)+ 

^ L (A.Fr.R[ra,Gn{i + l)-U{T b (i)-T a )\ V 

M.C p {-L(i + \).C p (T b (i)-l)-(u.A\(T b (i)-T.Ji + \))j 

With: 

The index ”i” specifies the variable's value at time i.At . We 
note from this equation that to find the temperature of the 
balloon at a given time, the conditions of weather at the 
considered instant and the temperature of the balloon at the 
moment iAt must be known, and other system temperatures 
are simultaneously calculated as follows: 
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r,(f+i)=r i (i+i)+(r,(i)-r t (i+i))e*p 


4 ^ 


R ,m.C pJ 


[4] 


T, (i + l)=r tat (f +1 )+{T, (i)-Tji + l))exp 


R t .m\C n 

\ t p J 


[5] 


T, o + 1)- T, 0 + 1) + 

rn .Cp 

[ 6 ] 


T(i + l)=Tjj +1)+ (7> (i)-T iBt (i + l))exp 


r L x ^ 
v R t m\C pJ 


[7] 


- The conductance between the balloon and the exterior 

of the absorber. 

- The resistance of the tubes between the base of the fin 

and the tube due to the solid fluid contact. 

- The length of the exchanger in the balloon and lengths 

of tubes between the balloon and the panel. 

- The flow rate of the heating fluid in the panel. 

- The flow rate of usage water. 

- The temperature surrounding the storage, the ambient 

temperature and the temperature of network. 

VII. RESULTS AND DISCUSSION 


V. THE GLOBAL SOLAR FLUX 

Any study or application of solar energy at a given site 
requires a complete knowledge and most detailed possible of 
the sunshine in site. This can be accomplished by the presence 
of a meteorological measuring station operating regularly for 
several years. In case of the absence of station, certain 
approximate methods for predicting characteristics of solar 
radiation are used. 

In our work a program in MATLAB to calculate the global 
solar radiation (direct, diffuse and reflected) is established 
Depending on weather conditions and pollution (climate, 
altitude, ...) at any point on the Earth at every instant; the 
latitude of the place, date (day of the year) and time in the day 
and at the angle of the wall exposure to the sun at the time 
concerned. 

The numerical code developed in MATLAB will be able to 
calculate the parameters involved using the values of solar 
radiation found in earlier calculations at a given location, for 
the best design of a solar heater namely influence of some 
parameters on the dynamic operation of the system by 
choosing the best settings according to the critical values of 
the efficiency of installation. 

VI. DATA 

A data file was selected from the database of our numerical 
code. The data come from daily weather measurements 
provides an annual file in hourly time. 

This file is relative to a given location: it consists of the 
following information: 

- The number of the day in each month when 

calculations are done. 

- The calculation step. 

- The Latitude of middle. 

- The orientation of the sensor relative to the south. 

- The solar flow. 

- The coefficient of albedo. 

- The proportion of insolation for 12 months. 

- The dimensions of the panel (length, width). 

- The thickness of the absorber. 

- The diameter of the tubes in the panel. 

- The coefficients of transmission, of absorption of the 

glazing. 

- The heat capacity of the heating fluid. 

- The dimensions of the storage balloon (radius, 

height). 

- The power of the pump. 


A. Influence of the orientation of the panel 

Generally, the orientation of the sensor must be in a position 
which provides a perpendicular direction of the incident solar 
irradiance 'Gn' on the panel plane. To determine the influence 
of the orientation of the panel, we varied the factor of this 
orientation "p", while other parameters are stabilized. For 
each p, we seek the efficiency of the installation rj. In figure 
(4), we showed the efficiency of the installation rj for different 
values of the angle of inclination of the panel. 

In our case, we considered that the case of Beirut of Latitude 
34,4°. This figure shows a maximum critical value for an 
efficiency orientation p ~ 42°. According to the general rule, 
the optimum orientation of the panel is equal to the latitude of 
+ 10°. Generally, in Beirut an optimum orientation of 42 °is 
considered. 


Figure 4: Evolution of efficiency for following the orientation 
of the panel. 



Figure 5 : Evolution of the storage temperature 
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B. . Influence of the ratio of the water mass of the balloon 
to the surface of panel (M/A) on the temperature of storage 
Tb 

A good choice of the ratio of the water mass of the balloon to 
the surface of panel (M/A) has a beneficial interest in 
improving the temperature of storage tank. 

In Figure (5), we show the evolution of the temperature of the 
storage during one day of the year for different values of the 
ratio [M/A]. We note that, from midnight, the temperature 
decreases due to heat loss through convection between the 
panel and the entourage then it increases with the sunrise and 
then it decreases after sunset. 

For the same surface of the panel, when the mass of water in 
the storage tank increases the same amount of heat must be 
transferred to a larger mass, therefore the storage temperature 
decreases. Thus it is beneficial to consider a larger collector 
surface for a better storage temperature. 

C. Influence of water mass ratio of the balloon on the 
panel surface 

The choice of a sensor's surface depends highly on the mass of 
water in the heat storage tank. The evolution of the efficiency 
of installation rj in function of ratio of water mass on the 
sensor surface is given in the figures (6 and 7). This evolution 
shows that for a given surface, when the mass of water 
increases, the sensor operates with a colder liquid thus it has a 
better efficiency till a critical value. Beyond this critical 
value, the efficiency improves only very few. In our case, the 
critical value (M/A) critical = 170 kg / m 2 . Using this critical 
value, for each panel's surface, the volume of the critical 
balloon can be obtained. 



Figure 6 : Evolution of efficiency in function of 


M/A (A = 10m 2 ) 



Figure 7: Evolution of efficiency in function of 
M/A (A = 6m 2 ) 


D. Influence of the mass flow of heat transfer fluid in the 

panel 

To determine this influence, the flow of heat transfer fluid in 
the sensor has been varied, while the other parameters are 
kept constant and then the efficiency of the installation rj is 
calculated. In Figure (8), the evolution of the efficiency of the 
installation in function of the coolant fluid flow is shown. For 
a given diameter, the amount of heat absorbed by the fluid 
increases with thus the efficiency increases. 

This evolution shows that the efficiency increases slightly for 
the values of flow less than 0,18 1/s then increases rapidly to a 
value 0,34 1/s and then increases slightly with the flow. 
Therefore, it is preferable to consider a flow rate of 0, 34 1/s. 

E. Influence of the heat capacity of the coolant 

As well to know the influence of the heat capacity of the 
coolant, this factor has been varied and then the efficiency of 
installation r\ is calculated. The figure (9) shows the evolution 
of efficiency of the installation for different values of the heat 
capacity of coolant. The amount of heat absorbed by the fluid 
increases with the heat capacity, therefore the efficiency 
increases .In addition, we notice that the water has the highest 
heat capacity (4180J / kg. ° C). This fluid is used despite the 
effects of corrosion and scaling which raises sustainable 
system operation problems. To solve the problems of 
corrosion effects, softeners or magnetic filters are used. 



Figure 8: Evolution of the efficiency of the installation 
according of mass flow rate of coolant fluid. 
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Figure 9: Evolution of the efficiency of the installation 
according to the heat capacity. 
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F. Influence of length of the heat exchanger in the storage 
balloon 

The heat exchange between the primary and the secondary 
circuits takes place through the heat exchanger in the storage 
tank. The improvement of the heat exchange is linked directly 
with the contact area, therefore the coil length L b should be 
precise to have a better efficiency of the installation. 
Obviously, this length depends on the volume of heat storage 
tank. 

As well to know the influence of the length of the exchanger in 
the storage tank, we fixed a volume for the balloon and we 
vary L b . The evolution of the efficiency of the installation in 
function to the latter brings us to a critical value of this length 
depending on the storage balloon volume. In the following, 
we give the coil length into the storage tank and the efficiency 
of the installation corresponding to five different values for 
the storage tank volume, which are: 150, 200, 250, 300, 350 
liters Figures (10, 11, 12, 13, and 14). 

The length of the exchanger in the storage tank increases as 
the energy transferred by the fluid to water is high, therefore 
the efficiency improves up until it reaches a critical value 
above which it becomes constant. In summarizing the five 
cases, a relationship can be concluded between the volume of 
the storage tank and the coil length: 
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Figure 10: V= 150 liters 
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Figure 11: V=200 liters 
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Figure 12: V= 250 liters 



Figure 13: V=300 liters 



Figure 14: V=350 liters 
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Figure 15: Variation of the length of the exchanger according 
to the volume of the balloon 
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Figure (15), shows the evolution of the length of the 
exchanger for different values of the volume of storage tank. 
Then from this figure, we can see an affine relation between 
the volume of the storage tank and the length of the exchanger 

by the following correlation: _ 

L b ~ 0,06 V + 21,4 
With L b in [m] and V in [m 3 ] 

G. Critical panel surface 

The increasing of the sensor's surface has a beneficial interest 
on the amount of heat absorbed, but this increase is 
accompanied with an increase in thermal losses with the 
environment. Thus the search for a critical surface is essential 
to improve the dynamic operation of the system as well as the 
economic aspect. To determine the critical surface of the 
panel, in the figure (16), the evolution of the efficiency of the 
installation depending on the panel surface is showed. It is 
clear that this efficiency increases greatly with the panel 
surface to a critical value. Beyond this critical value the 
efficiency decreases slightly, the decrease is due to the 
increased heat loss by convection between the sensor and the 
entourage. From this figure, we conclude that the critical 
surface is 10 m 2 . This critical value was found by Zelzouli et 
al [4]. 



Figure 16: Evolution of the efficiency of the installation 
depending on the panel surface. 


H. Energy requirements during the year 

We notice that the energy needed for heating and domestic hot 
water and during the first four months of the year and for the 
last two months, are high because the temperature and the 
solar flow are low during these periods, while during the 
intermediate month requirements are lower as shown in figure 
(17). Thus it is advantageous to use the extra energy to cover 
the lack of sunlight during the cold period. 



Year [month] 


Figure 17: Evolution of energy requirements.(V=250 m 3 ) 


VIII. CONCLUSION 

In this study, the most specific problems related to solar 
radiation and its conversion are treated. Our work enabled us, 
from different graphs to make a good design of a solar heater 
in a given location and this by selecting all critical parameters 
resulting in the highest efficiency or even when its variation is 
very little that the use of this parameter can save. 

We gave a theoretical analysis that can lead to tangible results 
following a work done on computer and essentially developed 
a numerical code able to calculate solar radiation in a given 
location, using the values found for calculating efficiency , 
storage power and other adopted parameters. 

To improve the dynamic operation of the system as well as the 
economic aspect, we investigated the influence of parameters 
on the efficiency of system. These results can be useful for 
designers of solar systems. 
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